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Aircraft design has been progressing toward reduced structural weight to improve 
fuel efficiency, increase performance, and reduce cost. Lightweight aircraft 
structures are more flexible than conventional designs and require new design 
considerations. Intelligent sensing allows for enhanced control and monitoring of 
aircraft, which enables increased structurally efficiency. The NASA Dryden Flight 
Research Center (DFRC) has developed an instrumentation system and analysis 
techniques that combine to make distributed structural measurements practical for 
lightweight vehicles. Dryden’s Fiber Optic Strain Sensing (FOSS) technology 
enables a multitude of lightweight, distributed surface strain measurements. The 
analysis techniques, referred to as the Displacement Transfer Functions (DTF) and 
Load Transfer Functions (LTF), use surface strain values to calculate structural 
deflections and operational loads. The combined system is useful for real-time 
monitoring of aeroelastic structures, along with many other applications. This paper 
describes how the capabilities of the measurement system were demonstrated using 
subscale test articles that represent simple aircraft structures. Empirical FOSS 
strain data were used within the DTF to calculate the displacement of the article and 
within the LTF to calculate bending moments due to loads acting on the article. The 
results of the tests, accuracy of the measurements, and a sensitivity analysis are 
presented. 

I. Introduction 

A. Value of Structural Shape and Load Measurement 

Knowledge of the operational shape and loading of a structure is valuable for a variety of aerospace 
applications, such as aeroelastic monitoring, improved control surface effectiveness, drag reduction, sonic 
boom reduction and stores pointing. 1 Existing camera-based methods are effective, 2 but often impractical 
for use as permanent installations for real-time data display on large, lightweight aircraft. 

The in-flight failure of Helios (Aerovironment, Inc., Monrovia, California), an extremely lightweight 
and flexible unmanned aircraft shown in Fig. 1, called attention to the need for a lightweight real-time 
deflection measurement system suitable for such an aircraft. Following the incident, the NASA Dryden 
Flight Research Center (Edwards, California) began investigating the use of distributed surface strain 
measurements as a means to calculate structural shape/ 1 

In addition to deflection measurement, the surface strain measurements of Fiber Optic Strain Sensing 
(FOSS) can be used to calculate the loads acting on a structure. 4 Current techniques use conventional 
metallic foil strain gages wired into shear, bending, and torsion bridges to measure loads at a few discrete 
points. 5 ' 6 Fiber Optic Strain Sensing and the Load Transfer Functions (LTF) enable load measurements 


1 Research Engineer, Aerostructures (RS) Branch, PO Box 273 M/S 48202A, and Member AIAA 

2 Research Engineer, Aerostructures (RS) Branch, PO Box 273 M/S 48202A, and Member AIAA 

1 

American Institute of Aeronautics and Astronautics 


distributed over the entire structure, often at a fraction of the weight of a conventional strain bridge 
arrangement. 



Figure 1. Helios (Aerovironment, Inc., Monrovia, California) aircraft; Note large operational wing 

deflections. 
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B. Equations 

A system of equations known as the Displacement Transfer Functions (DTF, also known as the Ko 
Displacement Theory) was developed at NASA Dryden. These equations use distributed surface strain 
values to calculate structural slopes and displacements. The DTF were formulated by integrating the 
curvature equation for the deformed-beam elastic curve. The structure is first discretized into multiple 
domains, so that the surface strain distribution within each domain can be described with a linear 
function. This piecewise linear approach enables integration of the beam curvature equation to yield slope 
and deflection equations in recursive formats. The slope and deflection equations were then combined 
into a single equation called the Displacement Transfer Function (Eq. 1). Furthermore, by calculating the 
deflection along a forward and aft fiber, the chord-wise angle of twist can be determined along the entire 
span. The formulation of the DTF is detailed in reference 7. 
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Similar to the DTF, a new method was developed at NASA Dryden to use distributed strain 
measurements to calculate operational loads using the LTF. This process requires the detennination of the 
cross-sectional properties of the structure by measuring strain values while a single known calibration 
load is applied. With these cross-sectional properties known at each domain, strain values induced by 
subsequent loads during operation can be used to determine the unknown loads acting on the structure. 
This paper considers loads only as bending moments; however, the method can be extended to include 
shear loads, as shown in Ref. 4. A flowchart of the LTF application process is shown in Fig. 2. Both the 
DTF and LTF equations are computationally efficient and are intended for use in a real-time environment. 
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Figure 2. Load transfer function implementation process. 
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C. Fiber Optic Strain Sensing Technology 

Fiber Optic Strain Sensing offers the ability to make thousands of distributed surface strain 
measurements with a lightweight, flight-worthy system (Fig. 3). Many different fiber optic 
instrumentation systems exist; however, the NASA-developed system is ideally suited to make distributed 
measurements on aerospace structures. The FOSS system used for the tests described in this paper was 
capable of making strain measurements at one -half-inch intervals along up to eight fibers, each up to 
40 feet in length at up to 100 samples per second. This allows nearly 8000 measurements covering 
320 linear feet. The FOSS systems under development further expand this capability. Fiber Optic Strain 
Sensing technology makes the DTF and LTF methods practical for flight vehicles. 



Figure 3. Fiber installation with strain gage shown for comparison. 

II. Test Description 


A. Test Articles 

The objective of the small-scale demonstration tests was to validate the shape and load measurement 
methods against empirical data on a variety of test articles under various loading conditions. The tests 
were designed to be simple, in order to investigate the underlying principles and limitations of the method 
using empirical data. Four 6061-T6 aluminum plates of different planform shapes served as the test 
articles (Fig. 4). Each of these plates was 0.190” thick with a span of 36 inches. The rectangular plate test 
article served as the baseline due to its simple geometry. The tapered plate and swept plate articles were 
each of shapes that are typical of aircraft wings. The open plate article was designed to have 
discontinuous cross-sectional properties, similar to a wing with internal ribs. One test article at a time was 
bolted into a cantilever fixture for testing, as shown in Fig. 4. 
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(c) Tapered Plate 


(e) Open Plate 


Figure 4. Test arrangement and test articles. 


B. Instrumentation 

Each of the test articles was instrumented with a single optical fiber which made several passes on 
both the upper and lower surfaces of the test article (shown in blue in Fig. 4). Strain measurements were 
made every one -half-inch along these fibers. 

The test articles were also instrumented with photogrammetry targets and a thennocouple, as can be 
seen in Fig. 5. The rectangular plate and swept plate also had metallic foil rosettes installed. The 
photogrammetry targets were placed along the fiber on the top surface of the test articles. An off-the-shelf 
two-camera photogrammetry system and its accompanying software were used during testing. The 
photogrammetric deflection measurements served as reference values to compare with the calculated 
deflections. The thermocouple provided a means to verify that the test article temperature did not change 
appreciably during testing. The strain rosettes provided a verification of the FOSS strain measurements. 
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Figure 5. Swept plate during uniform load case. 


C. Load Cases 

At least five unique load cases were applied to each test article. Each load case consisted of calibrated 
dead weights hung from different load points on the underside of the test articles, as shown in Fig. 6. The 
“uniform” load case was applied multiple times to each test article (called “Uniform A,” “Uniform B,” 
and “Uniform C”) to capture the repeatability of the test and calculations. The load cases applied to each 
test article are listed in Table 1. 
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Figure 6. Load points used during each load case are shown in solid red. 
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Table 1. Load cases performed on each test article. 



Uniform 

Tip- J 


load 

biased b 

Test article 

(A,B,C) 

load 


Rectangular plate 


X 


X 



Tapered plate 

X 

X 

X 

X 


X 

Swept plate 

X 

X 

X 

X 

X 

X 

Open plate 

X 

X 

X 

X 

X 

X 


Each load case was broken down into seven distinct load steps at which the test article motion was 
allowed to dampen out and photogrammetry images were captured. A description of the load steps is 
given in Fig. 7. 


Load Step 

Description 

1 

Unloaded 

2 

Minimum Load 

3 

Middle Load 

4 

Maximum Load 

5 

Middle Load 

6 

Minimum Load 

7 

Unloaded 



Figure 7. Numbering of load steps (hold points during testing). 


111. Data Processing 


A. Averaging at Load Steps 

Ten photogrammetry images were taken at each load step while the test article was in a static 
condition. The average of these ten images was used as the “true” vertical deflection for each target at 
each load step. The strain values were averaged over approximately 10 seconds at each load step. The 
effect of variation in the photogrammetric displacement values and strain values during these constant 
load conditions is discussed in the “Sensitivity Analysis” section. 

B. Interpolation and Extrapolation 

The spacing between the discrete domains used for deflection and load calculation was left to be 
variable, such that domains did not necessarily coincide with grating locations. Strain values at the 
domains were determined by linear interpolation using the two nearest measured strain values. 
Interpolation wasn’t possible at the root of the test article because it was inboard of the first grating; 
therefore, the strain values for the last two inches of each test article were extrapolated (as shown in Fig. 
8) using a fifth-order polynomial curve fit over the strain gratings along that fiber leg. For the purpose of 
this paper, the order of the polynomial curve fit was not formally optimized; a fifth-order curve fit was 
sufficient to characterize strains to the root. 
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Figure 8. Strain values at analysis domains were interpolated when between two adjacent strain 
measurements, and extrapolated to the root. 

C. Neutral Axis Calculation 

The method for determining the location of the neutral axis discussed in Ref. 7 was not used in the 
final data analysis. The referenced method uses the top and bottom surface strains to determine the 
distance from a fiber to the neutral axis, c, as shown in Eq. 2. Due to the quantity in the denominator, this 
method is unstable when the top and bottom strains are equal. 
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A new method was developed which eliminates the need to calculate the location of the neutral axis 
by changing the inputs to the DTF. This new method, developed by Dr. William Ko of NASA Dryden, 
uses the distance separating the top and bottom fibers, h, instead of the distance to the neutral axis, c, and 
uses the difference of the top and bottom strains, Sbottom- Stop, instead of strains from an individual surface, 
8 to p. The equivalency of this change is shown in Fig. 9. Thus, Eq. 1 becomes Eq. 3. This new method was 
found to eliminate errors associated with the instability in Eq. 2, and did not produce any adverse effects. 
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Figure 9. Change in strain location reference from neutral axis (a) to upper surface (b). 
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IV. Results 


A. Deflection 

The calculated deflection values generally compared very well with the photo grammetric 
measurements. Results for the middle fiber of the swept plate are shown in Figs. 10 and 11. Similar trends 
were observed in the other test articles and other load cases, which are shown in Appendix 1. Note that 
some load cases omit load steps 3 and 5 due to a lack of resolution in the available calibrated weights. 
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Figure 10. Swept plate deflection results for middle fiber during Uniform A load case. 
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Figure 11. Swept plate deflection results for middle fiber during single point load case. 


The deflection calculated along the trailing edge of the swept plate was consistently less accurate than 
in other areas. Toward the root of the swept plate test article, the principal strain direction quickly 
transitioned from the fiber direction to perpendicular to the root, as shown in Fig. 12. This is the area of 
greatest strain for the trailing edge fiber. The change in principal strain direction elsewhere is in areas 
with much less strain and did not cause the same inaccuracies. 
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Figure 12. Change in principal direction close to the root of the swept plate test article shown in 
finite element model. 

A summary of the tip deflection results is shown in Table 2. Comparing data at the maximum wingtip 
deflection (load step 4) from all test articles, all load cases, and all fiber legs, the root-mean-square error 
was 2.28 percent, or 0.106 inches. The greatest percent error was 10.9, which occurred along the trailing 
edge leg of the swept plate during the root-biased load case. The greatest absolute error was 0.312 inches, 
which occurred along the trailing edge leg of the swept plate during the trailing-edge -biased load case. 
Both of these maximum error cases are on the trailing edge leg of the swept plate test article, which as 
mentioned previously, was prone to inaccuracy. If this area is ignored, the greatest percent error was 4.7 
on the trailing edge leg of the open plate during the Uniform A load case, and the greatest absolute error 
was 0.229 inches on the trailing edge leg of the tapered plate during the leading edge biased load case. 
These accuracies (within 5 percent and .25 inches) were very encouraging for a strain-based method. 
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Table 2. Tip deflection calculation results. 



Rectangular Plate Deflection at Wingtip 


Trailing Edge 

Middle 

Leading Edge 


Calc 

Ref 

Ref - Calc 

Error 

Calc 

Ref 

Ref - Calc 

Error 

Calc 

Ref 

Ref - Calc 

Error 

Load Case 

(in) 

(in) 

(in) 

(X) 

(in) 

(in) 

(in) 

IX) 

(in) 

(in) 

(in) 

(X) 

Uniform A 

-5.488 

-5.398 

0.090 

-1.7 

-5.334 

-5.408 

-0.074 

1.4 

-5.360 

-5.408 

-0.048 

0.9 

Uniform 8 

-5.473 

-5.394 

0.079 

-1.5 

-5.309 

-5.404 

-0.095 

1.8 

-5.333 

-5.405 

-0.072 

1.3 

Uniform C 

-5.461 

-5.394 

0.067 

-1.2 

-5.212 

-5.403 

-0.191 

3.5 

-5.301 

-5.405 

-0.104 

1.9 

Single Point 

-2.202 

-2.172 

0.030 

-1.4 

-2.163 

-2.267 

-0.105 

4.6 

-2.282 

-2.358 

-0.076 

3.2 

Tip Biased 

-5.851 

-5.769 

0.081 

-1.4 

-5.727 

-5.774 

-0.047 

0.8 

-5.763 

-5.777 

-0.014 

0.2 

Root Biased 

-2.088 

-2.119 

-0.031 

1.5 

-1.962 

- 

- 

- 

-2.031 

-2.117 

-0.086 

4.1 

Leading Edge 

-4.014 

-4.035 

-0.021 

0.5 

-3.736 

-3.857 

-0.122 

3.2 

-3.730 

-3.676 

0.054 

-1.5 




Tapered Plate Deflection at Wingtip 


Trailing Edge 

Middle 

Leading Edge 


Calc 

Ref 

Ref - Calc 

Error 

Calc 

Ref 

Ref - Calc 

Error 

Calc 

Ref 

Ref - Calc 

Error 

Load Case 

(in) 

(in) 

(in) 

(X) 

(in) 

(in) 

(in) 

(X) 

(in) 

(in) 

(in) 

(X) 

Uniform A 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

Uniform B 

-6.493 

-6.507 

-0.014 

0.2 

-6.467 

-6.507 

-0.040 

0.6 

-6.379 

-6.519 

-0.141 

2.2 

Uniform C 

-6.448 

-6.512 

-0.064 

1.0 

-6.468 

-6.512 

-0.045 

0.7 

-6.376 

-6.524 

-0.148 

2.3 

Single Point 

-2.851 

-2.869 

-0.018 

0.6 

-2.863 

-2.881 

-0.018 

0.6 

-2.847 

-2.893 

-0.046 

1.6 

Tip Biased 

-7.110 

-7.141 

-0.031 

0.4 

-7.179 

-7.147 

0.032 

-0.4 

-7.121 

-7.161 

-0.040 

0.6 

Root Biased 

-2.355 

-2.396 

-0.041 

1.7 

-2.379 

-2.398 

-0.020 

0.8 

-2.351 

-2.403 

-0.052 

2.1 

Leading Edge 

-5.576 

-5.805 

-0.229 

4.0 

-5.777 

-5.763 

0.014 

-0.2 

-5.767 

-5.730 

0.037 

-0.6 



Swept Plate Deflection at Wingtip 


Trailing Edge 

Middle 

Leading Edge 


Calc 

Ref 

Ref - Calc 

Error 

Calc 

Ref 

Ref - Calc 

Error 

Calc 

Ref 

Ref - Calc 

Error 

Load Case 

(in) 

(in) 

(in) 

(X) 

(in) 

(in) 

(in) 

IX) 

(in) 

(in) 

(in) 

(X) 

Uniform A 

-7.361 

-7.661 

-0.300 

3.9 

-7.088 

-6.953 

0.134 

-1.9 

-6.229 

-6.233 

-0.005 

0.1 

Uniform B 

-7.361 

-7.655 

-0.294 

3.8 

-6.981 

-6.944 

0.037 

-0.5 

-6.172 

-6.219 

-0.046 

0.7 

Uniform C 

-7.358 

-7.658 

-0.300 

3.9 

-7.082 

-6.949 

0.132 

-1.9 

-6.244 

-6.225 

0.019 

-0.3 

Single Point 

-8.601 

-8.669 

-0.068 

0.8 

-7.876 

-7.695 

0.180 

-2.3 

-6.747 

-6.711 

0.036 

-0.5 

Tip Biased 

-6.409 

-6.576 

-0.167 

2.5 

-6.034 

-5.925 

0.110 

-1.9 

-5.259 

-5.255 

0.004 

-0.1 

Root Biased 

-2.059 

-2.310 

-0.251 

10.9 

-2.167 

-2.145 

0.022 

-1.0 

-1.919 

-1.976 

-0.057 

2.9 

Leading Edge 

-4.868 

-5.096 

-0.228 

4.5 

-4.745 

-4.714 

0.031 

-0.7 

-4.294 

-4.320 

-0.026 

0.6 

Trailing Edge 

-6.877 

-7.189 

-0.312 

4.3 

-6.481 

-6.447 

0.034 

-0.5 

-5.635 

-5.688 

-0.052 

0.9 




Open Plate Deflection at Wingtip 


Trailing Edge 

Middle 

Leading Edge 


Calc 

Ref 

Ref - Calc 

Error 

Calc 

Ref 

Ref - Calc 

Error 

Calc 

Ref 

Ref - Calc 

Error 

Load Case 

(in) 

(in) 

(in) 

IX) 

(in) 

(in) 

(in) 

(X) 

(in) 

(in) 

(in) 

(X) 

Uniform A 

-3.677 

-3.511 

0.166 

-4.7 

-3.499 

-3.513 

-0.015 

0.4 

-3.486 

-3.510 

-0.024 

0.7 

Uniform B 

-3.503 

-3.509 

-0.006 

0.2 

-3.497 

-3.513 

-0.015 

0.4 

-3.554 

-3.508 

0.047 

-1.3 

Uniform C 

-3.430 

-3.511 

-0.081 

2.3 

-3.483 

-3.514 

-0.031 

0.9 

-3.493 

-3.510 

-0.016 

0.5 

Single Point 

-3.430 

-3.401 

0.029 

-0.9 

-3.555 

- 

- 

- 

-3.738 

-3.722 

0.016 

-0.4 

Tip Biased 

-4.728 

-4.653 

0.075 

-1.6 

-4.719 

-4.649 

0.070 

-1.5 

-4.687 

-4.647 

0.040 

-0.9 

Root Biased 

-1.973 

-2.020 

-0.047 

2.3 

-2.002 

-2.020 

-0.017 

0.9 

-1.996 

-2.017 

-0.021 

1.0 

Leading Edge 

-4.430 

-4.595 

-0.165 

3.6 

-4.300 

-4.345 

-0.045 

1.0 

-4.062 

-4.092 

-0.029 

0.7 

Trailing Edge 

-4.171 

-4.096 

0.075 

-1.8 

-4.290 

-4.342 

-0.053 

1.2 

-4.574 

-4.592 

-0.018 

0.4 


Calc: Calculated value, using DTF and FOSS Strain 
Ref: Reference value, using photogrammetry data 


Ref-Calc: Difference between reference and calculated 
Error: (Ref-Calc)/Ref 

" - " : Value not available due to data corruption 
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B. Twist 

The angle of twist is calculated by comparing the deflection values at two chord-wise locations. 
Despite the accuracy of the individual deflection measurements, the calculated twist values were less 
accurate. The simple trigonometric relationship used to determine the angle of twist is very sensitive to 
errors in the calculated deflections. This caused errors of greater than 0.5 degrees at the tip of the test 
article in a majority of the load cases; results are shown in Table 3. The sensitivity of the twist equation is 
further examined in the “sensitivity analysis” section. 
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Table 3. Twist calculation results. 


Rectangular Plate Twist at Wingtip 

degrees) 

Load Case 

Calc 

Ref 

Ref- Calc 

Uniform A 

-0.67 

0.05 

0.72 

Uniform B 

-0.73 

0.06 

0.79 

Uniform C 

-0.83 

0.06 

0.89 

Single Point 

0.42 

0.97 

0.55 

Tip Biased 

-0.45 

0.04 

0.50 

Root Biased 

-0.30 

-0.01 

0.29 

Leading Edge 

-1.48 

-1.87 

-0.39 

Tapered Plate Twist at Wingtip (degrees) 

Load Case 

Calc 

Ref 

Ref- Calc 

Uniform A 

- 

- 

- 

Uniform B 

-2.17 

0.07 

2.24 

Uniform C 

-1.37 

0.06 

1.43 

Single Point 

-0.08 

0.13 

0.20 

Tip Biased 

0.21 

0.11 

-0.10 

Root Biased 

-0.07 

0.04 

0.11 

Leading Edge 

3.66 

-0.40 

-4.06 

Swept Plate Twist at Wingtip (degrees) 

Load Case 

Calc 

Ref 

Ref- Calc 

Uniform A 

-5.91 

-7.47 

-1.56 

Uniform B 

-6.20 

-7.51 

-1.31 

Uniform C 

-5.82 

-7.50 

-1.69 

Single Point 

-9.71 

-10.32 

-0.62 

Tip Biased 

-6.00 

-6.91 

-0.91 

Root Biased 

-0.73 

-1.74 

-1.01 

Leading Edge 

-2.99 

-4.04 

-1.04 

Trailing Edge 

-6.48 

-7.85 

-1.37 

Open Plate Twist at Wingtip (de 

grees) 

Load Case 

Calc 

Ref 

Ref - Calc 

Uniform A 

-0.99 

0.00 

0.99 

Uniform B 

0.27 

-0.01 

-0.28 

Uniform C 

0.33 

-0.01 

-0.34 

Single Point 

1.61 

1.68 

0.07 

Tip Biased 

-0.21 

-0.03 

0.18 

Root Biased 

0.12 

-0.02 

-0.14 

Leading Edge 

-1.91 

-2.62 

-0.71 

Trailing Edge 

2.10 

2.58 

0.49 


Calc: Calculated value, using DTF and FOSS Strain 
Ref: Reference value, using photogrammetry data 
Ref-Calc: Difference between reference and calculated 
" - " : Value not available due to data corruption 
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C. Load 

The load calculations produced promising results. The accuracy of the calculated loads, however, was 
found to be sensitive to the load distribution chosen for the calibration load, recalling that the calibration 
load is used to detennine the cross-sectional properties required for subsequent load calculations. During 
this study, one of the uniform load case applications served as the calibration load. Results from the load 
calculations closely followed the actual for subsequent uniform load applications. Calculations for 
different load cases follow the same trend as the actual values but were less accurate. Figures 13 and 14 
show the load calculation results for the tapered plate test article. Table 4 shows a summary of the load 
calculation results at the root of the test article during the maximum load step (load step 4) for each load 
case. The root-mean-squared bending moment error at the root is 10.7 percent; however, this number is 
greatly influenced by the large errors during the root-biased load cases. The error is much less for load 
cases that are similar to the calibration load, which indicates that using a realistic load distribution during 
the check loading will produce the best results during operation. 
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Figure 13. Tapered plate load calculation results for the middle fiber during the Uniform B load 
case. 
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Figure 14. Tapered plate load calculation results for the middle fiber during the leading edge load 
case. 
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Table 4. Summary of Load Transfer Function error; the first load case for each test article was 
used as the “calibration load.” 


Rectangular Plate Bending Moment at Root (in-lbs) 

Load Case 

Calc 

Ref 

Ref - Calc 

Error (%) 

Uniform A 

1200.0 

1200.0 

0.0 

0.0 

Uniform B 

1191.6 

1200.0 

8.4 

0.7 

Uniform C 

1156.1 

1200.0 

43.9 

3.7 

Single Point 

376.4 

385.0 

8.6 

2.2 

Tip Biased 

1034.2 

1176.0 

141.8 

12.1 

Root Biased 

720.5 

726.0 

5.5 

0.8 

Leading Edge 

835.6 

840.0 

4.4 

0.5 


Tapered Plate Bending Moment at Root (in-lbs) 

Load Case 

Calc 

Ref 

Ref - Calc 

Error (%) 

Uniform A 

- 

- 

- 

- 

Uniform B 

1200.0 

1200.0 

0.0 

0.0 

Uniform C 

1196.2 

1200.0 

3.8 

0.3 

Single Point 

403.1 

385.0 

-18.1 

-4.7 

Tip Biased 

1061.8 

1060.0 

-1.8 

-0.2 

Root Biased 

813.5 

726.0 

-87.5 

-12.1 

Leading Edge 

1145.8 

1110.0 

-35.8 

-3.2 


Swept Plate Bending Moment at Root (in-lbs) 

Load Case 

Calc 

Ref 

Ref - Calc 

Error (%) 

Uniform A 

480.0 

480.0 

0.0 

0.0 

Uniform B 

474.0 

480.0 

6.0 

1.2 

Uniform C 

474.9 

480.0 

5.1 

1.1 

Single Point 

374.5 

385.0 

10.5 

2.7 

Tip Biased 

327.4 

326.0 

-1.4 

-0.4 

Root Biased 

465.1 

352.0 

-113.1 

-32.1 

Leading Edge 

451.0 

408.0 

-43.0 

-10.5 

Trailing Edge 

393.0 

408.0 

15.0 

3.7 


Open Plate Bending Moment at Root (in-lbs) 

Load Case 

Calc 

Ref 

Ref - Calc 

Error (%) 

Uniform A 

480.0 

480.0 

0.0 

0.0 

Uniform B 

511.3 

480.0 

-31.3 

-6.5 

Uniform C 

450.8 

480.0 

29.2 

6.1 

Single Point 

361.0 

385.0 

24.0 

6.2 

Tip Biased 

498.0 

524.0 

26.0 

5.0 

Root Biased 

687.2 

522.0 

-165.2 

-31.6 

Leading Edge 

515.5 

600.0 

84.5 

14.1 

Trailing Edge 

597.8 

600.0 

2.2 

0.4 


Calc: Calculated value, using LTF and FOSS Strain 
Ref: Actual value 

Ref-Calc: Difference between reference and calculated 
Error: (Ref-Calc)/Ref 

" - " : Value not available due to data corruption 
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V. Sensitivity Analyses 


A. Deflection Sensitivity 

Following the initial deflection and load calculations, a sensitivity analysis was performed to 
determine the primary sources of error. The sources of error examined were strain values, plate thic kn ess 
values, and photogrammetry deflection values. The middle fiber of the swept plate test article during the 
Uniform C load case, load step 4, was used as the basis for the sensitivity study. 

The maximum observed standard deviation of a strain value at the investigated load step was 5.9 
microstrain. The deviation in strain values at several load steps was analyzed and determined to fit a 
normal distribution. For a normal distribution, the probability of the data falling within two standard 
deviations of the mean is approximately 95 percent. Twice the maximum standard deviation was applied 
to both the upper and lower surface strains (a total offset of 23.6 microstrain to the difference in strain 
values) as a conservative method for studying errors in the measured strain values. Errors in strain values 
may be attributed to several different causes. The most obvious cause for errors in strain values is 
measurement error in the Fiber Bragg Grating and FOSS interrogation system. Another cause for error in 
strain is incorrect location of the strain measurement. If the location of the strain measurement is 
incorrect, then the strain value is effectively incorrect. A final possible source for error is in poor fiber 
alignment during installation. Table 5 shows that errors in strain values can cause an error of just over 2 
percent for the swept plate. 

Table 5. Change in deflection calculation due to error in strain values. 



Offset 

Deflection 

Deflection 


(microstrain) 

(in) 

change (%) 

Measured 

n/a 

7.296 

n/a 

Measured + Offset 

23.6 

7.147 

-2.04 

Measured - Offset 

23.6 

7.445 

2.04 


Measured = lower surface strain - upper surface strain 


The thickness of the test articles was an input to the DTF and another potential source of error. For 
these tests, the nominal thickness (. 1 90 inches) of the aluminum plate was assumed to be constant. An 
offset of .015 inches was used for this analysis, which accounts for both the manufacturer’s thickness 
tolerance of the plate (±.009 inches) and the diameter of the fiber (.006 inches), which makes a significant 
difference in thickness for this test article. The deflection changes are nearly proportionally with 
thickness, so thickness errors are directly propagated into deflection errors, Table 6. 

Table 6. Change in calculated tip deflection due to change in thickness 
(performed using the swept plate test article, Uniform C load case, middle fiber). 



Thickness 

Thickness 

Deflection 

Deflection 


(in) 

error (%) 

(in) 

change (%) 

Nominal 

.190 

n/a 

7.296 

n/a 

Nominal + Offset 

.205 

7.89 

6.762 

-7.32 

Nominal - Offset 

.175 

-7.89 

7.921 

8.57 


The maximum standard deviation of a photogrammetry target displacement at the investigated load 
step was .0064 inches. An offset of twice this maximum standard deviation was applied to all the 
photogrammetry targets as a conservative method for studying the effect of errors in the photogrammetry 
values. This creates a range of ±.0128 inches around the measured value within which the true 
displacement can be expected to fall. Table 7 shows that photogrammetry uncertainty is insignificant (less 


19 

American Institute of Aeronautics and Astronautics 




than one percent) when measuring tip deflections. This uncertainty is assumed to be equal among all 
targets, which means that the photogrammetric error has the greatest influence at locations with small 
displacements. 

Table 7. Change in tip deflection due to error in photogrammetry values. 



Offset (in) 

Tip deflection 
(in) 

Tip deflection 
change (%) 

Measured 

n/a 

7.2511 

n/a 

Measured + Offset 

.0128 

7.2639 

0.18 

Measured - Offset 

.0128 

7.2383 

-0.18 


B. Twist Sensitivity 

The twist angle is calculated using the trigonometric relationship between the leading edge and 
trailing edge deflections, as shown in Eq. 4. This approach is very sensitive to errors in vertical deflection 
measurements (yLeading Edge and y Tr aiiing Edge)- Errors that would be considered low and acceptable for 
deflection can cause extremely large errors in the twist calculation. 


(j) = sin 


y LeadingEdg T Trai/ingEc&t 


* Chord 


( 4 ) 


Table 8 shows that 2.04 percent error in each of the deflection values (as calculated in the previous 
section) for the rectangular test article can cause an unacceptably large error (49 percent) in the calculated 
angle of twist. This extreme sensitivity to deflection errors makes the presented method for detennining 
twist angle impractical. 

Table 8. Change in twist calculation due to change in deflection. 


Reference Calculated Error 


Leading edge deflection 

Yi 

in 

6.25 

6.38 

2.04% 

Trailing edge deflection 

Yz 

in 

5.75 

5.63 

-2.04% 

Distance between fibers 

d 

in 

11.0 

11.0 

0% 

Calculated twist angle 

<J> 

deg 

2.61 

3.88 

49% 


C. Load Sensitivity 

The inputs to the LTF are strain values and loads measured during the calibration loading, and strain 
values measured during operation. The same variations in strain values described above were used in the 
analysis of the load equations (twice the maximum standard deviation). A potential variation in weight of 
.01 pounds was used. This corresponded to the resolution of the scale that was used to verify the 
calibration weights. Table 9 shows that these two sources of error could result in a measurement error of 
over six percent. 
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Table 9. Change in load measurement due to error in strain values and calibrated weights. 



Strain offset 
(microstrain) 

Load offset 
(lbs) 

Calculated 
moment (in-lbs) 

Change (%) 

Measured 

n/a 

n/a 

474.0 

n/a 

Measured + Offset 

11.6 

.01 

505.8 

6.7 

Measured - Offset 

11.6 

.01 

444.2 

-6.3 


VI. Discussion 


A. Applicability of Results 

The test articles used for this testing were simple, thin, monolithic and isotropic structures. Typical 
aerospace structures are much more complex. These tests, however, served as a valuable application of 
the measurement methods and provided an opportunity to collect and utilize empirical strain data within 
the DTF and LTF. The experiment uncovered limitations of the twist and load calculation methods, 
however the deflection calculation method was shown to be robust among the test articles and load cases. 

The DTF results from the described tests compare well with results from a test performed on a 
full-scale wing, 8 indicating that the method is not compromised by increasing test article size or 
complexity. 


VII. Conclusions 


The Deflection Transfer Functions were able to accurately measure wing deflection, with a 
root-mean-squared error of 2.28 percent of the actual deflection. These results were seen as a success for 
the technology. The DTF equations are not particularly sensitive to moderate errors in strain values. The 
relatively thin geometry of the plate test articles used for this testing made the equations more sensitive to 
errors in thickness measurement. The 7.89-percent thickness error used in the sensitivity analysis would 
correspond to nearly one-half inch on a 6-inch-thick wing; in reality, wing thickness can be measured 
much more accurately. 

The twist measurement method proved to be poorly suited to use with empirical data due to a great 
sensitivity to error in the deflection measurements. This sensitivity makes the twist measurement method 
impractical unless if the error in the deflection measurement is greatly reduced. 

Though not yet as accurate as photography-based deflection measurements, these strain-based 
methods offer unique capabilities. For many applications, the benefits of the weight reduction and real- 
time capability will outweigh the reduced accuracy. 

The LTF results are also very encouraging. The root-mean-squared error value listed above (10.7 
percent) is greater than a conventional strain bridge system may exhibit (in certain circumstances lower 
than 3 percent error); however, many of the load cases exhibited less error than a conventional system. 
The load measurement method would be particularly valuable when the load distribution that will be 
experienced in operation is understood; a calibration load should be chosen which simulates the load 
distribution that the structure will encounter during operation. 
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APPENDIX 1: Deflection Plots 
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Figure Al l. Rectangular plate deflection results for leading edge fiber during Uniform A load case. 
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Figure Al_2. Rectangular plate deflection results for middle fiber during Uniform A load case. 
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Rectangular Plate -Uniform A Load Case -TE Fiber 


A * m * » » , 

* A ‘ A # 




=§. -2 

c 

O 

o 

<D 

S - 3 


-4 


-5 


-6 


• Step 1 , Calculated 

• Step 2, Calculated 

• Step 3, Calculated 

• Step 4, Calculated 

Step 5, Calculated 

x- Step 6, Calculated 

Step 7, Calculated 
Step 1 , Photogrammetry 
Step 2, Photogrammetry 
Step 3, Photogrammetry 
Step 4, Photogrammetry 
Step 5, Photogrammetry 
Step 6, Photogrammetry 
Step 7, Photogrammetry 


■ ■ * 


10 


15 20 

Span (in) 


25 


30 


35 


Figure Al_3. Rectangular plate deflection results for trailing edge fiber during Uniform A load case 
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Figure Al_4. Rectangular plate deflection results for leading edge fiber during Uniform B load case 
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Rectangular Plate -Uniform B Load Case -Mid Fiber 
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Figure Al_5. Rectangular plate deflection results for middle fiber during Uniform B load case. 


Rectangular Plate -Uniform B Load Case -TE Fiber 


l t 1 




*- » 


E -2 

c 

O 


a - 3 


-4 


-5 


-6 


» ft 


V 


Step 1 , Calculated 
Step 2, Calculated 
Step 3, Calculated 
Step 4, Calculated 
Step 5, Calculated 
Step 6, Calculated 
Step 7, Calculated 
Step 1 , Photogrammetry 
Step 2, Photogrammetry 
Step 3, Photogrammetry 
Step 4, Photogrammetry 
Step 5, Photogrammetry 
Step 6, Photogrammetry 
Step 7, Photogrammetry 


10 


15 20 

Span (in) 


25 


30 


35 


Figure Al_6. Rectangular plate deflection results for trailing edge fiber during Uniform B load case. 
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Rectangular Plate -Uniform C Load Case -LE Fiber 
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Figure Al_7. Rectangular plate deflection results for leading edge fiber during Uniform C load case. 


Rectangular Plate -Uniform C Load Case -Mid Fiber 


* * I s ; * * 


E -2 

C 

O 


a - 3 


-4 


-5 


-6 


A 


A 


V 


Step 1 , Calculated 
Step 2, Calculated 
Step 3, Calculated 
Step 4, Calculated 
Step 5, Calculated 
Step 6, Calculated 
Step 7, Calculated 
Step 1 , Photogrammetry 
Step 2, Photogrammetry 
Step 3, Photogrammetry 
Step 4, Photogrammetry 
Step 5, Photogrammetry 
Step 6, Photogrammetry 
Step 7, Photogrammetry 


* 




10 


15 20 

Span (in) 


25 


30 


35 


Figure Al_8. Rectangular plate deflection results for middle fiber during Uniform C load case. 
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Figure Al_9. Rectangular plate deflection results for trailing edge fiber during Uniform C load case. 
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Figure A110. Rectangular plate deflection results for leading edge fiber during single point load case. 
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Figure Al ll. Rectangular plate deflection results for middle fiber during single point load case. 
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Figure Al_12. Rectangular plate deflection results for trailing edge fiber during single point load case. 
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Rectangular Plate -Tip Biased Load Case -LE Fiber 
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Figure A113. Rectangular plate deflection results for leading edge fiber during tip-biased load case. 
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Figure A114. Rectangular plate deflection results for middle fiber during tip-biased load case. 
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Figure A115. Rectangular plate deflection results for trailing edge fiber during tip-biased load case. 
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Figure A116. Rectangular plate deflection results for leading edge fiber during root-biased load case. 
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Figure A117. Rectangular plate deflection results for middle fiber during root-biased load case. 
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Figure Al_18. Rectangular plate deflection results for trailing edge fiber during root-biased load case. 
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Figure A119. Rectangular plate deflection results for leading edge fiber during leading edge load case. 
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Figure Al_20. Rectangular plate deflection results for middle fiber during leading edge load case. 
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Figure A121. Rectangular plate deflection results for trailing edge fiber during leading edge load case. 
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Figure A122. Tapered plate deflection results for leading edge fiber during Uniform B load case. 
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Figure A123. Tapered plate deflection results for middle fiber during Uniform B load case. 


Tapered Plate -Uniform B Load Case -Trailing Edge Fiber 


* ft ft v# it ft ft * *- » A i * ■ * jfc: a a A A 

^ A ft vfc ^ ft • . 

. K • R ft a 


-2 


f -3 


CD 

Q 


-5 


-7 


• Step 1, Calculated 

• Step 2, Calculated 

• Step 3, Calculated 

• Step 4, Calculated 

Step 5, Calculated 
Step 6, Calculated 
Step 7, Calculated 
Step 1 , Photogrammetry 
Step 2, Photogrammetry 
Step 3, Photogrammetry 
Step 4, Photogrammetry 
Step 5, Photogrammetry 
Step 6, Photogrammetry 
Step 7, Photogrammetry 


• p ft ». & ». ». 


ft ft _ 


10 


15 20 

Span (in) 


25 


30 


35 


Figure Al_24. Tapered plate deflection results for trailing edge fiber during Uniform B load case. 
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Figure Al_25. Tapered plate deflection results for leading edge fiber during Uniform C load case. 
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Figure A126. Tapered plate deflection results for middle fiber during Uniform C load case. 
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Figure A127. Tapered plate deflection results for trailing edge fiber during Uniform C load case. 
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Figure Al_28. Tapered plate deflection results for leading edge fiber during single point load case. 


36 

American Institute of Aeronautics and Astronautics 


-0.5 


-1.5 


-2 


-2.5 


Tapered Plate -Single Point Load Case -Mid Fiber 

• • • • • •. • • • • V . 

• Z * • 


A 


A 

A 


V 


Step 1 , 
Step 2, 
Step 3, 
Step 4, 
Step 5, 
Step 6, 
Step 7, 
Step 1 , 
Step 2, 
Step 3, 
Step 4, 
Step 5, 
Step 6, 
Step 7, 


Calculated 

Calculated 

Calculated 

Calculated 

Calculated 

Calculated 

Calculated 

Photogrammetry 

Photogrammetry 

Photogrammetry 

Photogrammetry 

Photogrammetry 

Photogrammetry 

Photogrammetry 


A 


10 


15 20 

Span (in) 


25 


30 


35 


Figure Al_29. Tapered plate deflection results for middle fiber during single point load case. 
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Figure A130. Tapered plate deflection results for trailing edge fiber during single point load case. 
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Figure Al_31. Tapered plate deflection results for leading edge fiber during tip-biased load case. 
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Figure A132. Tapered plate deflection results for middle fiber during tip-biased load case. 
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Figure Al_33. Tapered plate deflection results for trailing edge fiber during tip-biased load case. 
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Figure A134. Tapered plate deflection results for leading edge fiber during root-biased load case. 
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Figure Al_35. Tapered plate deflection results for middle fiber during root-biased load case. 
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Figure A136. Tapered plate deflection results for trailing edge fiber during root-biased load case. 
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Figure A137. Tapered plate deflection results for leading edge fiber during leading edge load case. 
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Figure A138. Tapered plate deflection results for middle fiber during leading edge load case. 
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Figure A139. Tapered plate deflection results for trailing edge fiber during leading edge load case. 
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Figure A140. Swept plate deflection results for leading edge fiber during Uniform A load case. 
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Figure A141. Swept plate deflection results for middle fiber during Uniform A load case. 
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Figure A1 


42. Swept plate deflection results for trailing edge fiber during Uniform A load case. 
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Figure A143. Swept plate deflection results for leading edge fiber during Uniform B load case. 
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Figure A144. Swept plate deflection results for middle fiber during Uniform B load case. 
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Figure A145. Swept plate deflection results for trailing edge fiber during Uniform B load case. 
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Figure A146. Swept plate deflection results for leading edge fiber during Uniform C load case. 
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Figure A147. Swept plate deflection results for middle fiber during Uniform C load case. 
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Figure A148. Swept plate deflection results for trailing edge fiber during Uniform C load case. 
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Figure A149. Swept plate deflection results for leading edge fiber during single point load case. 
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Figure A150. Swept plate deflection results for middle fiber during single point load case. 
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Figure Al_51. Swept plate deflection results for trailing edge fiber during single point load case. 
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Figure A152. Swept plate deflection results for leading edge fiber during tip-biased load case. 
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Figure A153. Swept plate deflection results for middle fiber during tip-biased load case. 
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Figure A154. Swept plate deflection results for trailing edge fiber during tip-biased load case. 
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Figure A155. Swept plate deflection results for leading edge fiber during root-biased load case. 
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Figure A156. Swept plate deflection results for middle fiber during root-biased load case. 
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Figure Al_57. Swept plate deflection results for trailing edge fiber during root-biased load case. 
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Figure A158. Swept plate deflection results for leading edge fiber during leading edge load case. 
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Figure A159. Swept plate deflection results for middle fiber during leading edge load case. 
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Figure A160. Swept plate deflection results for trailing edge fiber during leading edge load case. 
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Figure A161. Swept plate deflection results for leading edge fiber during trailing edge load case. 
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Figure A162. Swept plate deflection results for middle fiber during trailing edge load case. 
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Figure Al_63. Swept plate deflection results for trailing edge fiber during trailing edge load case. 
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Figure A164. Open plate deflection results for leading edge fiber during Uniform A load case. 
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Figure A165. Open plate deflection results for middle fiber during Uniform A load case. 
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Figure A166. Open plate deflection results for trailing edge fiber during Uniform A load case. 
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Figure A167. Open plate deflection results for leading edge fiber during Uniform B load case. 
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Figure A168. Open plate deflection results for middle fiber during Uniform B load case. 
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Figure A169. Open plate deflection results for trailing edge fiber during Uniform B load case. 
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Figure A170. Open plate deflection results for leading edge fiber during Uniform C load case. 
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Figure A171. Open plate deflection results for middle fiber during Uniform C load case. 
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Figure A172. Open plate deflection results for trailing edge fiber during Uniform C load case. 
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Figure A173. Open plate deflection results for leading edge fiber during single point load case. 
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Figure A174. Open plate deflection results for middle fiber during single point load case. 
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Figure A175. Open plate deflection results for trailing edge fiber during single point load case. 
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Figure Al_76. Open plate deflection results for leading edge fiber during tip-biased load case. 
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Figure A177. Open plate deflection results for middle fiber during tip-biased load case. 
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Figure A178. Open plate deflection results for trailing edge fiber during tip-biased load case. 
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Figure A179. Open plate deflection results for leading edge fiber during root-biased load case. 


Open Plate -Root Biased Load Case -Mid Fiber 



L L 

L L L L ■ : 


*1 { ; 1 ■ * * Cf r 

* ft- >* : * • , . 


’ * * * 

? ft • • ? # 


* . ft » 

*■ * R • • sj* 


• Step 1 , Calculated 

* A . 


• Step 2, Calculated 

• 9 ' • 


• Step 3, Calculated 



• Step 4, Calculated 

• ^ “ • 

• -- ft • 


Step 5, Calculated 

• ^ ft • 

■* ’■ • # 


Step 6, Calculated 

* ^ % • 

• . ,u. 


Step 7, Calculated 

• 

• 


Step 1 , Photogrammetry 

• __ 

• 


Step 2, Photogrammetry 

• 

/ft 


Step 3, Photogrammetry 

• 


Step 4, Photogrammetry 

• 


Step 5, Photogrammetry 

• 


Step 6, Photogrammetry 



Step 7, Photogrammetry 

r r r r 


0 5 10 15 20 25 30 35 

Span (in) 


Figure A180. Open plate deflection results for middle fiber during root-biased load case. 
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Figure A181. Open plate deflection results for trailing edge fiber during root-biased load case. 
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Figure A182. Open plate deflection results for leading edge fiber during leading edge load case. 
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Figure A183. Open plate deflection results for middle fiber during leading edge load case. 
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Figure A184. Open plate deflection results for trailing edge fiber during leading edge load case. 
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Figure A185. Open plate deflection results for leading edge fiber during trailing edge load case. 
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Figure A186. Open plate deflection results for middle fiber during trailing edge load case. 
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Figure Al_87. Open plate deflection results for trailing edge fiber during trailing edge load case. 
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Figure A2_l. Rectangular plate load results for middle fiber during Uniform A load case. 
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Figure A2_2. Rectangular plate load results for middle fiber during Uniform B load case. 
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Figure A2_3. Rectangular plate load results for middle fiber during Uniform C load case. 


400 


350 


300 


250 


200 


150 


100 


Rectangular Plate - Single Point Load Case - Mid Fiber 


-50 


Step 1 , Calculated 
Step 2, Calculated 
Step 3, Calculated 
Step 4, Calculated 
Step 5, Calculated 
Step 6, Calculated 
Step 7, Calculated 
Step 1 and 7, Actual 
Step 2 and 6, Actual 
Step 3 and 5, Actual 
Step 4, Actual 



10 


15 20 

Span (in) 


25 


30 


35 


Figure A2_4. Rectangular plate load results for middle fiber during single point load case. 
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Figure A2_5. Rectangular plate load results for middle fiber during tip-biased load case. 
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Figure A2_6. Rectangular plate load results for middle fiber during root-biased load case. 
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Figure A2_7. Rectangular plate load results for middle fiber during leading edge load case. 
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Figure A2_8. Tapered plate load results for middle fiber during Uniform B load case. 
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Figure A2_9. Tapered plate load results for middle fiber during Uniform C load case. 
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Figure A210. Tapered plate load results for middle fiber during single point load case. 
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Figure A211. Tapered plate load results for middle fiber during tip-biased load case. 
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Figure A212. Tapered plate load results for middle fiber during root-biased load case. 
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Figure A2_13. Tapered plate load results for middle fiber during leading edge load case. 
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Figure A214. Swept plate load results for middle fiber during Uniform A load case. 
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Figure A215. Swept plate load results for middle fiber during Uniform B load case. 
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Figure A216. Swept plate load results for middle fiber during Uniform C load case. 
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Figure A217. Swept plate load results for middle fiber during single point load case. 


Swept Plate - Tip Biased Load Case - Mid Fiber 


350 


300 


250 


200 


150 


100 


-50 


Step 1, Calculated 
Step 2, Calculated 
Step 3, Calculated 
Step 4, Calculated 
Step 5, Calculated 
Step 6, Calculated 
Step 7, Calculated 
Step 1 and 7, Actual 
Step 2 and 6, Actual 
Step 3 and 5, Actual 
Step 4, Actual 



10 


15 20 

Span (in) 


25 


30 


35 


Figure A218. Swept plate load results for middle fiber during tip-biased load case. 
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Figure A219. Swept plate load results for middle fiber during root-biased load case. 
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Figure A2 20. Swept plate load results for middle fiber during leading edge load case. 
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Figure A2 21. Swept plate load results for middle fiber during trailing edge load case. 
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Figure A2 22. Open plate load results for middle fiber during Uniform A load case. 
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Figure A2 23. Open plate load results for middle fiber during Uniform B load case. 
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Figure A2 24. Open plate load results for middle fiber during Uniform C load case. 
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Figure A2 25. Open plate load results for middle fiber during single point load case. 
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Figure A2_26. Open plate load results for middle fiber during tip-biased load case. 
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Figure A2 27. Open plate load results for middle fiber during root-biased load case. 


600 


500 


400 


300'*- 


200 


Open Plate - Leading Edge Load Case - Mid Fiber 


Step 1, Calculated 
Step 2, Calculated 
Step 3, Calculated 
Step 4, Calculated 
Step 5, Calculated 
Step 6, Calculated 
Step 7, Calculated 
Step 1 and 7, Actual 
Step 2 and 6, Actual 
Step 3 and 5, Actual 
Step 4, Actual 


100 


-100 



10 


15 20 

Span (in) 


25 


30 


35 


Figure A2_28. Open plate load results for middle fiber during leading edge load case. 
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Figure A2 29. Open plate load results for middle fiber during trailing edge load case. 
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